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Abstract 

The string equation for the [P, Q}= Q formulation of non-perturbatively stable 2D 
quantum gravity coupled to the (2m — 1,2) models is studied. Global KdV flows 
between the appropriate solutions are considered as deformations of two compatible 
linear problems. It is demonstrated that the necessary conditions for such flows 
to exist are satisfied. A numerical study reveals such flows between the pole-free 
solutions of pure gravity (m = 2), the Lee- Yang edge model (m = 3) and topological 
gravity (m = 1). We conjecture that this is the case for all of the m-critical models. 
As the m = 1 solution is unique these global flows define a unique solution for each 
m-critical model. 



April 1992. 



Introduction 

It is well knownfl] [2] that the non-perturbative sector of the [P, Q] = 1 
formulation of the (2m — 1,2) minimal models coupled to gravity is problematic 
for m even. The string equations for the string susceptibility p of this definition are 
the Painleve I hierarchy [1], (m+ \)t m TZ m [p\ = z. For m even, the solutions for p 
with the asymptotics p — * z x l m for z — > ±oo are the complex truncated Boutroux- 
type solutions[3], which are unacceptable. (Without loss of generality we may 
set t m = l/(a m (m+ ^)), where a m is the coefficient of p m in the Gel'fand-Dikii 
polynomial 1Z m .) Real solutions with the z — > +oo asymptotic above have poles, 
which are also unacceptable [2]. An analysis of the Dyson-gas picture of the matrix- 
model realisation of these critical models [4] reveals a fundamental instability in the 
critical matrix eigenvalue configurations for m even. 

The models have an underlying KdV-flow structure, dt m p = 7Z m+1 , which 
describes the evolution of p towards the mth model in the interpolating string 
equation Dm=i( m + ^)t m T^m = z. Many attempts were made to surmount the 
problem of the non-perturbative instabilty of pure gravity (m = 2) as defined by the 
Painleve I equation. One of these [5] was to try to define the solution for the string 
susceptibility by flowing from the stable m = 3 solution[6], switching on the massive 
coupling £2- The result of this numerical work showed the flow to be unstable at 
large t^'- The solution for the string susceptibility develops wild oscillations in the 
negative z region, signalling the tunnelling of eigenvalues in the matrix model and 
presumably heralding the onset of the generic singular solution of Painleve I in 
the limit. In ref.[7] the theory of isomonodromic deformations was used to study 
such flows between solutions of the string equations with the required asymptotics. 
This study showed analytically that there is no path between the m = 2 and m = 3 
models via the KdV flows. This is true for all neighbouring m-critical models in 
this definition[8]. The KdV flows, although present infinitesimally in the Painleve I 
definition, fail to materialize globally. 

Later, another definition of the above (2m — 1,2) models was proposed, 
which shared the perturbative results of the previous definition, but lacked the 
non-perturbative instabilities, as analysis of its original matrix model definition 
shows[9][10]. This is the [P, Q] = Q definition, the central principle of which is the 
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non-perturbative preservation of the KdV flows. This results in the most general 
string equations compatible with the KdV flows which are 

uR 2 - -nn + \{nf = o (i) 

for the string susceptiblity u, where TZ = J2m=i( m + |)^m^-m — z for the inter- 
polating case and we set t m = 2/(a m (2m + l))5 m ,o for the pure m-critical model. 
The solution for u for such a model has the asymptotics u — > z x l m ', for z — > ±oo. 
These boundary conditions can be shown to fix the integration constants so that 
the number of solutions at each m-critical point is discrete [10]. Numerical studies 
have revealed [10] [11] just one such solution for each of the m = 1, 2 and 3 models 1 . 
In each case the solution was manifestly pole-free. 

These results suggest the existence of a unique, pole-free solution to each of 
the m-critical models, thus completing the program of presenting a definition of the 
(2m — 1,2) models coupled to gravity which is stable everywhere. In this letter we 
wish to address this suggestion. There exist powerful techniques 2 for the analysis 
of equations such as (1). This is due to the fact that they are related to non-linear 
evolution equations which are solvable via the Inverse Scattering Transform[15], 
in this case the KdV heirarchy. Indeed, as discussed in ref.[16], the equation (1) 
for the mth critical model defines (a sub-class of) the self-similar solutions of the 
mth member of the KdV heirarchy. Recall that the string equations found for the 
unitary matrix models [17] are the self-similar solutions of the mKdV hierarchy. The 
non-perturbative solutions for the family of critical models for that hierarchy have 
been shown to exist [18] [19]. This was carried out by the use of a Gel'fand-Levitan- 
Marchenko-type linear integral equation [20] [21] equivalent to the string equation. 
Standard inverse scattering techniques are then used, resulting in the definition of 
a one-parameter family of pole-free solutions for each model, where for one value 
of the parameter, there is a solution with the required asymptotics. 

We would like to carry out a similar analysis for our string equations (1). We 
derive an integral equation equivalent to (1), but due to technical difficulties to 

1 In the case of m = 1, a simple change of variables[10] maps to the problem of 
proving that there exists a unique pole-free solution of the Painleve II equation. 
Such a solution is known to exist [12] [13]. 

2 See for example ref.[14]. Other references will appear later in this letter. 
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be described later it is not clear that using it we can construct arguments which 
facilitate a proof of uniqueness and freedom from poles. Instead, in this letter we 
take a different route. Following the lead of refs.[14] [7], we formulate our string 
equation as a flatness condition for two first order linear systems. The KdV flows 
in this framework can be shown to represent isomonodromic deformations of one 
of these linear systems and isospectral deformations of the other. We study the 
form of the monodromy data when we impose the asymptotics required for our 
solutions for the string susceptiblity. We find that there are no global obstructions 
to KdV flows between the m-critical models. Turning to a numerical study, we 
demonstrate that such KdV flows are indeed present by tracking the flows between 
the m = 1, 2 and 3 solutions. This then allows us to define a unique solution for the 
string susceptibility of each m-critical model by global KdV flow from the unique 
m = 1 solution. 



The analysis of the string equation. 

The starting point is to rewrite the differential of (1) as the compatibility 
condition [L,P] =0 for two linear problems: 

— A r). 

d fo ^ 



PV = P=-— + 

dX V 7 -« 

where 7 = j^(z + J2T=o ^ k ^m- k)i a = 7 /2, and (3 = 7(A + it) — a The compatibility 
condition above arises naturally in the context of isomonodromic deformations, 
reviewed in[22]-[23]. We follow the methods in those works closely, and refer the 
reader to them for details. The method concentrates on the second equation in (2). 
This equation is singular at the origin and at infinity and we diagonalize the most 
singular part at infinity, setting 

V^W, A=(*\ -_\?) (3) 

We make the change of variable £ 2 = A, transforming the equation into 

dW 
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MW (4) 



where now ( are the Pauli matrices) M = (l/2£ — 2a()ai + i{(3 — / yC 2 ) <J 2 + (P + 
7C 2 )c"3. An asymptotic treatment around C = 00 reveals that there exists a formal 
solution ^ of (4) [22] [24] 

$ = (1 __^_ (ri + ...)e(< 2m+ V 2 ( 2m + 1 )+^)-3 (5) 

where the dots represent higher order terms and diagonal parts, which are increas- 
ingly (with m) complicated functions of u and its derivatives. The formal solution 
is such that if we divide up a neighbourhood of 00 into sectorial domains fi^, then 
in each domain there is a true solution of (4), ^k, which is asymptotic to (5). 
Furthermore, on O/c+i nfijt, ^fc+i = ^kSk for some triangular matrix Sk called 
a Stokes matrix. (We will also use Sk to label the overlap flk+i Hfifc.) In our 
case there are 4m + 2 sectors around infinity defined by (k = 1, 2, ... , 4m + 2) 
fi fc = [7r/2 + 7t(/c-2)/(2m + 1) < arg( < tt/2 + nk/{2m + 1)]. Each sector S k 
contains a unique ray, 7r(/c + m)/2m + 1, called a Stokes line along which one of 
the fundamental solutions is maximally dominant. The initial lines of the sectors 
Sk are called conjugate Stokes lines, denoted Ck- The Stokes matrices are generally 
triangular [22] and, taking into account the sign of the leading term 



>2fc+l 



( 1 0\ _(l s 2k 

U*+i J' S2k ~\o 1 



As our equation is singular at the origin we should in principle do a similar 
analysis of the solution in a neighbourhood of this singularity. However, the original 
singularity before the change of variables is regular, so the analysis is not so 
complicated. We return to the original equation (2) and make the substitution 
A = C 2 and because the singularity at the origin is not diagonalizable we split the 
solution matrix $ into its components $>\ and i = 1, 2 and study the equation 
in terms of these. Keeping only highest order terms in l/( we easily find that 
the fundamental solutions are $i j2 = (b((), a(() + b(()ln() where a(£)and b(() are 
holomorphic in a neighbourhood of the origin. The monodromy property of this 
solution round the origin is then 

*(Ce 2 ™) = $(C)J, J=(l 2 l l ). (6) 



In order to compare the solutions at different points we need to introduce 
connection matrices relating $ and These are defined by \& = where \1/ = 
and A was defined in equation (3). The monodromy data parametrizing the 
solutions are thus the Stokes matrices Sj, the monodromy matrix round the origin, 
J, and the connection matrix, C. Not all of these parameters are independent, and 
the number of independent parameters may be reduced by studying the symmetries 
of the matrix M, and imposing reality on the string susceptibility u. Further 
relations are obtained by studying the behaviour of the solutions as we encircle the 
origin once. Half of the 4m + 2 original Stokes parameters are determined by the 
symmetries and one by the behaviour around the origin, leaving 2m parameters. Of 
the four parameters in the connection matrix two are determined by the behaviour 
around the origin, one by the requirement detC = 1 and the last one describes a 
global scaling freedom in the solution [14]. 

In ref.[7] the corresponding problem for the Painleve I hierarchy was studied 
resulting in necessary conditions for the global KdV flows to exist between m = 2 
and m = 3 and these conditions were proven not to be satisfied. In short the 
requirement is that the sectors in which the Stokes matrices are non-trivial should 
overlap between the different models, because the KdV How equations preserve the 
monodromy data[7] [14] . The Stokes matrices can be calculated by the method in 
ref.[23]. However we will not do so here as there is a pictorial way to proceed in 
order to derive the necessary conditions. In ref.[7] it was proven that when four 
conjugate Stokes lines emerge from a point, the Stokes matrix corresponding to 
the middle sector must be trivial, i.e. unity. So we will impose our boundary 
conditions and use this argument to find the sectors that must be trivial. First 
we derive necessary conditions on the Stokes matrices so that u ~ z L l m as z — > oo. 
Keeping only highest order terms as z — > oo we rescale, change variables £ = At 1 / 2 ™, 
u = r 1 / 7 ™ and determine the components of equation (2) as r — > oo: a = 0, (3 = r(l + 
A 2 )1/2A 2 £ X 2k a m - k , and 7 = r 1 " 1 / m l/2A 2 £ \ 2k a m - k . The a m are the coefficients 
of u m in 1Z m . As in the appendix of ref.[23] we write the equation in component 
form 

d\ ua 

where P = r/r(A) and Q = (1 + A 2 ) £ £ a k a j X 2( - 2m ~ k ~ j ~ 1 ' ) . Here r(A) is also a 
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polynomial but its explicit form is of no interest to us and VP 1 ' 2 are the vector 
components of This equation can be solved by standard WKB methods and the 
approximate solution is 

*1&kb) = (P 2 Q)- 1/4 exp(±r 1 + 1 /2 T n| Q i/2^ ) 

The important part here, for determining necessary conditions, is simply the form 
of the polynomial Q. It always has a double-zero at the origin, a simple zero at 
X = ±i and 21 — 2 distinct zeros of order two elsewhere. We may determine the 
lay of conjugate Stokes lines by taking into account their behaviour at infinity 
and the fact that they only meet at zeros of / Q x l 2 dz (they are defined through 
ReJQ 1 / 2 dz = 0). We can now use the pictorial argument of [7] to determine that 
the Stokes matrices corresponding to sectors Si with i odd and 1 < i < 2m + 1 are 
trivial 3 . 

If we now consider solutions which decay rapidly enough as z — > — oo the matrix 
M can be replaced by M = (x + l/2( 2m )as and the solutions are simply ^l/ 1 ' 2 = 
exp(±(£ 2m+1 /2(2m + 1) + zQ). Using this, we determine the conjugate Stokes 
lines and by repeating the argument of ref.[7] we find that the Stokes matrices 
corresponding to sectors Si with i even and 1 < i < 2m + 1 must be trivial 3 . In 
this way we find that if the string susceptibility u is to behave asymptotically as 
desired then it is necessary that the nontrivial monodromy is concentrated along 
the imaginary axis. The only non-trivial Stokes matrices are Si, SWre+i, 5 < 2m+2 and 
5 < 4m+2- that this is only a necessary condition. For example for the m=l equation in 
the unitary matrix case the coresponding boundary conditions actually completely 
determine the Stokes parameters [23]. In figs.l, 2 and 3 we display the configurations 
of conjugate Stokes lines for the cases m = 1, 2 and 3. 

Thus the necessary condition for the existence of flows between the models, 
derived in ref.[7], is satisfied here. We will later confirm their existence by the 
numerical study of the flows between the m = 1, 2 and 3 models. 

Note that the WKB solutions above and the structure of the Stokes sectors are 
exactly the same as for the Painleve II hierarchy [18]. In fact it will be shown in a 

3 This result becomes a conjecture for m > 3, as in ref.[7] due to the increasing 
complexity of the configuration of Stokes lines. 
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later work [16] that the existence of a Muira map between the two heirarchies allows 
further study of the solutions to equation (1). 



The Integral Equation. 

We now consider briefly the inverse problem to that above, i.e. given the Stokes 
matrices we can compute u, deriving an integral equation which can be studied by 
the methods of refs.[12] [21]. 

We immediately specialise to solutions with the asymptotics that we desire. 
Thus, as derived before we have only the Stokes matrices Si, S^m+i, 5 < 2m+2 and 
(the sectors closest to the imaginary axis) non-trivial. An important 
observation is that because of the singularity at the origin the monodromy around 
this point is not trivial (as it is for the Painleve II heirarchy) but is given by 
equation (6). A consequence of this is a branch cut along the positive imaginary 
axis. Therefore ^ *i+4 m +2 (superscripts denote components and subscripts 
sectors). Another consequence is that the relations among the Stokes parameters 
which follow from the symmetries of the equation are more complicated than in the 
Painleve II case. The relations are s\ + S4 m +2 = s = —2i and Sk = Sk+2m+i for the 
Stokes parameters The main idea 4 is to use Cauchy's theorem and relate an 
integral of ty 1 along the path C\ to one along C<i (see fig. 4) and to continue this 
around the singular point £ = oo. Thus we consider the integral 



Jo 



where 9 = —(( 2m+1 /2(2m + 1) + x() and for simplicity we assume that the pole ( 
lies in sector Si. We proceed by using Cauchy's theorem 

The first term on the right hand side is the residue at £ = £, the second term arises 
from an integral at infinity and 71 is a segment of a circle around the origin. Next we 



JC^ t-C 



4 See, for example. ref.[14] for details. 
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use the relations among the solutions in neighbouring sectors Vf^ = ^ ~ s i^ < 2^ 
to rewrite the third term in equation (8) 



ic 2 c-c ->c 2 f-c 

Note that ^^ e ® ~* 00 as C ~~ > 00 i n ^2 and therefore we can not relate the term 

* (2) e^ 

— Si J Ca d£ by contour integration to an integral along C3. This integral 

appears in this form in the final equation. We can however express the other 
integral along C2 in terms of an integral along C3 

1 fr* = -3 ( ) + k tr« + L tr« (9) 

where now ^ = vl/jj . We then continue in this manner around infinity and sum 
all the resulting equations to find 

*( 1 )(C)e^O = 

1\ 1 f V^e 9 8l f ^e 6 „ s f ^e 6 „ ^ 



OJ 2niJc(-£ 2ni Jc a C - £ 2ttz Jc b (-£ 

where = &*\ The main difference between this equation and the corresponding 
equation in case of the Painleve II hierarchy is that the integrals along C\ and C4 m+3 
(on the two sides of the branch cut along the imaginary axis) no longer cancel. The 
contours C a (Cb) run inward along C2 m +2 then clockwise (anticlockwise) around the 
origin and out along C2 (C^m+s). The contour C starts at £ = 00 travels to the left 
of the branch cut, encircles the origin and returns to infinity along the other side of 
the branch cut. A corresponding equation can be derived for the other component 

tfW(C)e-*(C) = 

f°^__L f y (2)e ~ 9 g i f ^ (1) e~ g s Sl r ^ f(2) e~ 6> ( U ) 
V 1 J ~ 2ri Jc (-£ ~ 2ri Jc c (-£ ^ 2riJc b (-£ 

where C c runs inward along C2 m +3 and out along C^ m ^. Thus equations (10) 
and (11) are the integral equations which correspond to our string equation for 
the solutions with the asymptotics u(z) ~ 2 1 / 7n ,0 as z — > ±00. (By comparison 
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with ref.[14] we find that these equations are precisely the same as in the case of 
the Painleve II heirarchy with a non-zero constant on the right hand side. This 
is a further consequence of the existence of a map between that heirarchy and the 
string equations (1). The significance of these results will be discussed in ref.[16].) 
In principle these integral equations can now be analysed using the methods in 
[14] [12] [13]. Due to the complications from the singularity at the origin, we find it 
difficult to proceed via this route. To study the uniqueness of the required solutions, 
we turn once again to the study of the KdV flows. 

The existence of global KdV flows. 

The flow between solutions of two critical points m and n is governed by 
the massive interpolating string equation which is (1) with TZ = (m + \)t m lZ m + 
(n + \)t n 1Z n — z. In flowing from m to n, t m is redundant and we may rescale 
variables to set it to l/(a m (m + |)), where a m is the coefficient of u m in 7Z m . 
t = (n + \)a n t n defines the one-parameter flow d t u ~ 7l' n+1 . Starting at the 
solution u m (z) for the pure m-critical theory we evolve t from zero, and track 
the solution u(z, t) to the interpolating equation, using the fact that the solution 
can only change infinitessimally via tha KdV flows [10]. Scaling this solution to 
u(z,t) = t~ 2 / 2n+1 u(zt~ l / 2n+1 ,t), we arrive at the solution to the pure n-critical 
model in the limit t — > +oo: u n (z) = limt^ +00 u(z, t). 

This is precisely the procedure that we carried out. We set up the interpolating 
string equation as a two-point boundary value problem for arbitrary t. We solved 
the first differential of the string equation as this allows for better numerical 
behaviour. Otherwise it is necessary to calculate 1/TZ to evaluate the highest 
derivative in (1), which results in rounding errors for large positive z. We ensured 
that we obtained solutions to the correct equation (1) by comparison with the 
pure m-critical solutions obtained in previous numerical treatments of this string 
equation 5 . Here, we calculate the positive z boundary conditions at arbitrary t by 
solving the tree-level equation u m + tu n = z numerically. The negative z boundary 
condition was is u = throughout. 

5 The pure m = 2 solution is displayed in ref.[10], and the pure m = 3 in ref.fll]. 
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We set up the flow from m = 3 to m = 2 and from m = 2 to m = 1. At each 
t, the differential equation was solved after typically 12 iterations, using the NAG 
FORTRAN library routine D02RAF, with a global error estimate of ~ 10~ 7 in the 
case ofm = 2^m = l and ~ 10~ 4 for m = 3 — > m = 2. Some of the intermediate 
stages are displayed in figs. 5 and 6, starting at t = and going up to t = 150. For 
large t (>100) the routine required an increasingly higher number of iterations to 
meet its error tolerances. However upon comparison of u(z, t) as defined above 
with the pure solutions it was observed that the limit procedure described above 
was converging rapidly; well before any numerical problems due to very large t could 
develop. 

The above convergence of the procedure is in striking contrast to the flow in 
ref.[5] using the 7Z m = z definition. There the instabilty of the flow at large t 
showed that the KdV flows could not connect the solutions in the limit. Here, we 
have successfully demonstrated the opposite: Global KdV flows exist between the 
pole-free m = 1, 2 and 3 models. We find it a reasonable conjecture that this is true 
for all the m-critical models. 

The existence of these global flows defines for us a unique solution to each of 
the m-critical models by flow from the unique m = 1 solution. 
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Figure Captions 



Fig. 1: The m = 1 Stokes sectors for (a) u ~ z as z — > +00 and (b) it ~ as 



-00 



Fig. 2: The m = 2 Stokes sectors for (a) tt~z2 as z — > +00 and (b) u ~ as 



-00 



Fig. 3: The m = 3 Stokes sectors for (a) « ~ z3 as z — > +00 and (b) u ~ as 

2 — >• —00 

Fig. 4: The integration contours used to derive the integral equation 

Fig. 5: Some snapshots of the evolution from the m = 2 (pure gravity) model (top 
curve) to the m = 1 model. The values of the KdV-flow parameter are 
0.0, 0.5, 1.0, 2.0, 5.0, 10.0, 30.0, 50.0, 70.0 and 150.0. A rescaling of the 
final curve, as described in the text, shows rapid convergence to the pure 
m = 1 solution. 

Fig. 6: Some snapshots of the evolution from the m = 3 (critical Lee- Yang) model 
(top curve) to the m = 2 model. The values of the KdV-flow parameter 
are 0.0, 0.5, 1.0, 2.0, 5.0, 10.0, 30.0, 50.0, 70.0 and 150.0. A rescaling of 
the final curve, as described in the text, shows rapid convergence to the 
pure m = 2 solution. 
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